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Abstract 

The density of two initially independent condensates which are allowed to expand and overlap 
can show interferences as a function of time due to interparticle interaction. Two situations are 
separately discussed and compared: (1) all atoms are identical and (2) each condensate consists of 
a different kind of atoms. Illustrative examples are presented. 
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The investigation of interferences between particles is one of the most basic tools to learn 
on the nature of quantum gases. Interferences attracted much attention in particular in the 
case of Bose- Einstein condensates (BECs) both from the theoretical and experimental sides, 
see, e.g., ) la In a popular set up studied, identical atoms are produced in two 

traps which we may call the left and right traps and which are separated by a barrier. By 
removing the traps and the barrier between them, the atoms expand freely and can overlap. 
In experiment, the photographs obtained show spectacular interference fringes [l, 2]. 

The interference of two parts of a single coherent condensate is by now well understood, 
see, e.g., 4], si, Q]. On the other hand, relatively little is known on the interference of two 
initially independent (i.e., fragmented) BECs, except for the case of non-interacting particles 



Fragmented BECs can be produced using a barrier between the two traps which is 
so high and broad that tunneling between them is negligible. 

In the available experiments, the atoms are prepared in a double-well trap potential and 
it is not generally proven whether the atoms form a coherent BEC, a fragmented BEC, or 
a combination thereof. However, it is feasible nowadays to produce in the lab two spatially 
separated, initially independent BECs, see, e.g., [10(, and this allows for experiments with 
definitely fragmented BECs. Apart from its importance as a fundamental problem, the 
solution of the problem of interference of two initially independent condensates is thus also 
of practical relevance. 

In the scenario of two initially independent BECs the initial state of the many-body 
system before removing the traps reads 

\^) = (N L \N R \)- l/2 ^ L ) NL ^ R f R \vac), N L + N R = N, (1) 

where the b\ and b R are the usual creation operators for bosons in the left and right traps, 
respectively, which contain definite numbers and N R of atoms in them. After removing 
the traps, the state \^f) is no longer an eigenstate of the system's Hamiltonian Ho and 
expands in space as a function of time. The time-dependent density, i.e., the expectation 
value of the density operator p(x) as a function of time becomes [?J 

p{x, t) = (*(t) \p{x)\ *(t)> = N L \$ L {x, t)\ 2 + N R \$ R (x, t)\ 2 , (2) 

where the $t t R(x, t) are the single-atom states corresponding to b L)R {t) = 
exp(iH t)bL tR exp(—iH t). Obviously, the density is a sum of the individual densities of 
the two condensates and does not exhibit an interference term. 
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We would like to draw attention to the fact that the literature result (J2D has been obtained 
under the assumption that atoms belonging to the two different BECs do not interact with 
each other. Very recently it has been demonstrated that, in the presence of interaction, the 



density p(x, t) does show an interference term 
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12|: 



p(x, t) = Pll(x, t) + p RR (x, t) + Plr(x, t), (3) 

where pll and prr are the densities of the expanding separated BECs as if the two BECs 
do not communicate, and plr is the change of the density due to the interaction between 
them. The terms contributing to Pll (prr) contain only &l(&r) and 6^(6^) operators, e.g., 
and those contributing to plr contain only mixed products, e.g., &£&] ? &l&r. The 
finding ([3]) has many consequences. In particular, the corresponding interference structures 
remain after the statistical averaging over many experimental runs. Of course, as p(x, t) 
changes in time, the average must be carried out at the same value of t. 

To derive ([3]) the full Hamiltonian H = Hq + V of the system after removing the traps 
including the particle-particle interaction V has been taken into account. For the ease of 
presentation, we employ the widely used contact interaction V(x,x') = X5(x — x'), where 



A is proportional to the s-wave scattering length 0, 



l^ . Of course, any other interparticle 



interaction can be used as well. As usual, H describes the motion of the free atoms. Starting 
from H = Hq + V and the initial many-body state (CQ), we have obtained the exact result 
for p(x, t) up to first order in the particle-particle interaction strength A. The corresponding 
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expression is somewhat lengthy and is not given here, but can be found in 
briefly mention properties of this result. Clearly, the interference term piR vanishes for 
t — *■ 0. Furthermore, piR{x,t) vanishes as expected if the atoms do not interact with each 
other (A — ► 0). The interference term piR{x,t) is enhanced by the product N^Nr of the 
numbers of atoms in the two initial BECs. 

The above discussions make clear that the interaction between the particles gives rise to 
an interference term in the density of two initially independent BECs of identical bosons. 
Before presenting a numerical example we go one step further and pose the question whether 
we can formulate a mean-field theory which reproduces exactly the many-body small A 
result mentioned above. Such a theory would open the door for real applications. The 
usual mean-field theory leads to the well-known and widely used Gross-Pitaevskii equation 
which reproduces exactly the density of BECs in a coherent state in the weak interaction 
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limit [7|, |l_3j]. Clearly, this equation is inapplicable to BECs in fragmented states ((Tj). For 
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gmented states a more general multi-orbital mean-field theory has been recently derived 
14] . In the present scenario two orbitals are involved and the respective time-dependent 
mean-field [TDMF(2)] takes on the appearance (for the general derivation of TDMF, see 



iil> L = V 



h + \(N L - 1) |^ L | 2 + 2XN R \?p R \ 2 \ ip L , 

h + X(N R - 1) |^| 2 + 2AiV L \^ L \ 2 } ip R (4) 



where the initial conditions are ipL,R{x,t = 0). h is the usual one-particle Hamiltonian (in 
our scenario just the kinetic energy operator — ^J^) and V — 1 — \ipL )(^t\ — \ip R ){^r\ is 
a projector which ensures orthonormalization of the orbitals ipi and ip R [151 ] . In TDMF(2) 
the density can be expressed by p(x,t) = Nl t)\ 2 + N R \i/)r(x, t)\ 2 . Indeed, it can 

be shown [3] that the TDMF (2) exactly reproduces the many-body result in the weak 
interaction limit. 

In the following we apply the TDMF(2) theory (j3J). For coherent states the time- 
dependent Gross-Pitaevskii equation, which is exact in the weak interaction limit, has been 
demonstrated in many cases to be applicable for intermediate and stronger interactions 

nn 

[3, H3|. Similarly, there is reason to expect that for fragmented states the TDMF theory, 



which has been proven to be exact in the weak interaction limit |l2j , is applicable well beyond 
this limit. We mention that TDMF(l) is nothing but the time-dependent Gross-Pitaevskii 
equation. 

We consider harmonic traps centered at ±Xo, each containing a BEC with interaction 
A = 0.1. At t — these traps are removed. As initial conditions iPl,r{x,0) we choose the 
respective solutions of the stationary Gross-Pitaevskii equation at this A to account for the 
interaction when the harmonic traps are released. In Fig. [1] the density p(x, t) computed 
using the TDMF(2) equations is shown as a function of time. As seen in the figure, at t = 
the density consists of two separated distributions centered at ±xo- The traps are removed 
at this time and the distributions start to broaden and to overlap. At about t = 3 one begins 
to see impact of the interference term in the density which becomes strongly pronounced as 
time proceeds. 

We see that the density of two initially independent condensates which are allowed to 
overlap can show interference effects in the presence of interparticle interaction. The physics 



of so called fragmented states, like the state in Eq. (pQ), is generally very different from that of 



coherent states 16j . Coherent states of condensates have been extensively studied, mostly in 
the framework of the Gross-Pitaevskii equation [?], [3] . A BEC in a coherent state can exhibit 



interference fringes even in the absence of interaction [7, [9), 0, [3]. Take, for instance, 
the coherent state \^ coh ) = (Nl)- 1 / 2 ^)* \vac) with W = (b[ + b ] R )/^2. This immediately 
leads to p coh (x, t) = y |$r,(x, t) + t)| 2 in the absence of interaction between the atoms, 

and hence to the interference term pfj^ix, t) = NKe ($2^/?)- F° r expanding Gaussians with 
initial width 2a located at ±Xo, the oscillatory part of p c £ R is simply given by cos[K(t)x] with 
K(t) = 8xo(t/m)/(a 4: + 4t 2 /m 2 ). This interference term is qualitatively different from that 
arising due to the interaction between the particles. Another important difference between 
p c £ft(x,t) and pLn{x,t) worth mentioning is that the former depends on the relative phase 
between <3>£ and while the latter does not depend on this phase. 

Whether in an experiment the initial state is coherent or fragmented depends on the 
experimental conditions. It is beyond the scope of this work to argue whether or not the 
initial state in the currently available experiments on interference is fragmented. It is also 
not our intention to take side in the ongoing debate on whether these experiments detect 
the density or higher-order correlation functions, although we tend to share the opinion of 
some researchers see, e.g., , Il7| . that the density is measured. What we can state, 



is that if one measures the density of two freely expanding initially independent BECs, 
it will only show interferences in the presence of interaction. This leads to the following 



the fact that nowadays one can vary the 
Two measurements are necessary. 



proposal for an experiment which makes use o 
strength of the interaction between the atoms [19 
If the measurement with interaction shows interferences which disappear upon measuring 
with the interaction turned off, then (a) the initial state was a fragmented state and (b) the 
interaction is responsible for the interferences. 

Until now the indistinguishability of the atoms has been considered a precondition for 
interference effects. In experiments on interferences one often starts with a single coherent 
BEC made of identical bosons and produces two BECs by ramping up a barrier. Interferences 
are then observed after removing the traps and the barrier. Nearly all theoretical works on 
interferences in BECs rely on the property of coherence of the whole system consisting of 



identical particles. Even the two very recent works 
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12 1 which discuss interferences of 



identical independent BECs due to interaction between the particles have assumed that 
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indistinguishability is necessary. We show below that this is unnecessary. 

As above, we consider two initially independent BECs in an initial state like in ([1]). Now, 
however, each of these two BECs is made of a different kind of atoms. For simplicity we 
call them "left" and "right" atoms, and assign to them the creation operators b L and b R in 
(PQ). As usual, the Hamiltonian H = H + V now contains three interactions terms V = 
Vl + Vr + Vlr accounting for the interaction between the "left" particles [Vl = XlS(x — x')], 
between the "right" particles [Vr = Xr6(x — x')], and between the particles of both kinds 
[Vlr = XlrS(x — x')}, respectively. Using the same basic techniques as in [12j for identical 
particles, we obtain equation (jHJ) for the density also in the present case of distinguishable 
BECs. The density contains a term plr due to the interaction Vlr between the particles of 
the two BECs. Moreover, the analytic expression for p in the weak interaction limit is very 



similar to that in 



12] for identical particles. 



In contrast to the case of indistinguishable bosons, where the time-dependent mean-field 
equations for fragmented condensates have been derived very recently Jl5l|, time-dependent 



mean-field equations for mixtures of different bosons are well known 21 
present situation the latter read 



22 



23|. For the 



h + X L (N L - 1) |^| 2 + XlrNr \^ r \ 2 
h R + Xr(N r - 1) |^| 2 + X LR N L |Vl| 2 



IpR 



(5) 



where for simplicity we have used the same nomenclature for the orbitals ipi,R{x,t) as for 
identical particles and the kinetic energies Jil,r = — 2mL - differ due to the possibly 
different masses of the left and right bosons. It is not surprising that (jSJ) reproduce the 
exact result for the density in the weak interaction limit. 

It is illuminating to briefly compare equations (jSJ), which we call time-dependent coupled 
Gross-Pitaevskii [TDCGP(2)] equations, to the TDMF(2) for identical bosons (0}. For 
identical bosons there are only a single mass m and interaction strength A, and a factor 
2 appears in (j3j) due to the exchange of identical bosons. More importantly, TDMF(2) 
maintains the orthogonality of the orbitals ipL and ipR while TDCGP(2) does not. 

To demonstrate that the density exhibits an oscillatory pattern also for interacting dis- 
tinguishable condensates, we show a few examples computed via (jSJ). For simplicity we put 
m L = m R = m. In our first example we choose Xlr = 2Xl = 2Xr = 0.2 which leads to the 
analogous scenario discussed in Fig. [1] for identical atoms. At t — the initial density is as 



in Fig. [TJ At later times the density can evolve differently than in Fig. [TJ only because the 
orbitals ipL arid ipR do not have to be orthogonal to each other for distinguishable atoms. 
In Fig. [2] the density is shown for t = 8 and compared with the analogous density obtained 
with (J3J) for identical particles. Both identical and distinguishable systems exhibit oscilla- 
tory structure but their differences are substantial. In the lower part of Fig. [2] we show the 
individual subdensities Nl \?Pl\ 2 and Nr \iPr\ 2 at t = 8 and remark that even a moderate 
overlap (ipL \ iPr) can have a considerable impact on the oscillatory pattern of the density. 

In our next example we choose Xl = Xr = and Xlr = 0.2, implying that the atoms 
in each of the left and right condensates do not interact, but those belonging to different 
condensates do. At t = we thus have two normalized Gaussians localized at the minima 
±Xo of the harmonic traps. After removing the traps these Gaussians expand and overlap, 
and an oscillatory pattern develops. The result is shown for t = 22 in the upper panel 
of Fig. [3j Until now the interactions studied were repulsive and we address the question 
whether an oscillatory structure can also arise for attractive interactions. To answer this 
question we investigate the same scenario but with Am = —0.2. The result is depicted in the 
middle panel of Fig. [3j Remarkably, the oscillatory structure is even much more pronounced 
than for the repulsive interaction. 

In our last example we ask whether an oscillatory pattern always evolves when interactions 
are present. We now choose Xl = Xr = Xlr = 0.1 and find no oscillations in the density up 
to quite long times. The density at t — 10 is shown in the lower panel of Fig. [3] This result 
can be understood from (jSJ). Assuming uil = uir, Xl = Xr = Xlr and Nl = Nr 3> 1, 
we notice that both equations (jSJ) are actually equivalent as they depend only on the total 
density p = Nl\i^l\ 2 + N R \iPr\ 2 . The analogous situation is not possible for identical 
particles because of the factor 2 appearing in (j3j) due to the exchange interaction. 

The origin of the oscillations in the density is interference. Commonly, one attributes a 
phase difference to the appearance of interferences. This is particularly simple in the case of 
two parts of a coherent condensate which show interference effects when overlapping. One 
may attribute a phase to each of these parts. As discussed above for fragmented condensates 
of identical particles, the relative phase of the fragments is irrelevant in the context of 
interferences, and this is, of course, also the case for condensates made of different kinds of 
particles. In these situations, interparticle interactions are responsible for interference effects. 
Here, we may speak of interaction-assisted self-interference. Consider for simplicity two 
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freely expanding, initially non-overlapping, different condensates with repulsive interaction 
between them. In each of these condensates the orbital has a phase which depends on x 
and t. As usual for a freely expanding isolated condensate, its orbital and hence its phase 
are smoothly changing such that no interference occurs. Once the interacting condensates 
begin to overlap, atoms are decelerated in the overlapping region and the local phase changes 
there. At later times the changed part of the orbital is superposed with other expanding 
parts of the same coherent condensate and this leads to interferences. In the extreme case of 
infinitely strong repulsive interaction between the two condensates, each of the condensates 
is reflected from the other condensate as if it were an expanding hard wall. The interferences 
then arise from the superposition of the reflected and advancing parts of the orbital. 

The theory presented here is easily extendable to any kind of interparticle interaction. It 
is also easily extendable to the case where one does not let the two BECs expand freely by 
removing the traps completely. One may, e.g., remove only the barrier and let the BECs 
expand in the new global trap. Since the interference structures depend on the interaction, 
a wealth of effects can be expected by varying the interaction, the form of the individual 
traps and the numbers Nl and Nr of the particles. In the case of distinguishable particles 
also the difference in masses rriL, vtir and, most importantly, the intra- and inter-condensate 
interactions \l, \r and \lr enrich the possible range of interference phenomena. 
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FIG. 1: (Color online) The density p(x,t) of two condensates of 500 atoms each for A = 0.1 as a 
function of time computed with TDMF(2) (black) compared to the density pli + prr of two BECs 
which do not interact with each other, each computed with the Gross-Pitaevskii equation (red). 
The quantities shown are dimensionless. For more details see text. 



10 




-40 -30 -20 -10 10 20 30 40 

X 



FIG. 2: (Color online) The density of two condensates each made of 500 atoms of a different kind 
computed with TDCGP(2) for Xl = = 0.1 and Xlr = 0.2 at t = 8 (black). For comparison 
also the analogous result for identical particles computed with TDMF(2) (for details see Fig. [T]) is 
shown (red). The lower panel shows the respective subdensities Nl \ iPl\ 2 for distinguishable atoms 
(black) and Nr \iPr\ 2 for identical particles (red). 
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FIG. 3: The density of two condensates each made of 500 atoms of a different kind. Upper panel: 
repulsive interaction of the two condensates, Xl = Xr = 0, Xlr = +0.2. Middle panel: attractive 
interaction of the two condensates, Xl = Xr = 0, Xlr = —0.2. Lower panel: Xl = Xr = Xlr = 0.1. 
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